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Résumé 

This article aims to make explicit the characteristic projection introduced by 
Panchishkin in [7j and gives various results on how to calculate it efTectively on 
a computer. It ends with concrète examples putting those ideas in application. 

1 Introduction 

On souhaite dans cet article montrer comment rendre explicite et effective infor- 
matiquement la méthode de la projection caractéristique introduite par Pantchichkine 
dans [7], qui sert à obtenir des congruences pour les valeurs spéciales des fonctions L 
de formes modulaires. Pour cela, on va dans cette introduction commencer par fixer 
quelques notations sur la théorie des formes modulaires, et rappeler les énoncés prin- 
cipaux nécessaires. On expliquera ensuite ce qu'est la projection caractéristique, puis 
on énoncera les résultats proprement dits. 

La suite de l'article consistera en des preuves et explications plus précises et quelques 
exemples concrets de calculs (avec l'aide de SAGE0). 

Pour les définitions précises et les preuves, le lecteur est invité à consulter des 
références comme les livres de Serre [8] , de Miyake [6] ou de Diamond et Shurman [4] . 
On trouvera aussi un survol plus détaillé et tourné vers les calculs effectifs dans le livre 
de Stein [9] : c'est en effet là que sont expliqués les algorithmes utilisés par SAGE pour 
calculer les espaces de formes modulaires et les opérateurs de Hecke associés. 

Etant donnés k ^ 2 et N ^ 1, A4k(Ti(N)) dénote l'espace vectoriel des formes 
modulaires de poids k et de niveau N pour Ti(N), et Sk(T 1 (N)) son sous-espace 
des formes paraboliques, sur lequel on va se concentrer ; on peut voir un élément / 
soit comme une fonction z i— > f(z) avec z dans le demi-plan de Poincaré, soit via le 
développement en série de Fourier à l'infini comme une fonction q \-ï f(q) avec q dans 
le disque unité, ce second point de vue étant privilégié. 



1. version 5.4, disponible sur http//www . sagemath. org 
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On sait en outre que ces espaces peuvent se décomposer en somme directe de sous- 
espaces, indexée par les caractères de Dirichlet \ de conducteur divisant N : 

f MkpiiN)) = ® x M k {N lX ) 
\ SkÇT^N)) = ® x S k (N, x ) 

Entre ces espaces de formes modulaires, on sait définir des familles d'opérateurs ; en 
particulier les opérateurs de Hecke, (T„)„> 2 . La famille est engendrée par les opérateurs 
(T p ) p premier, qui commutent deux à deux. Lorsque p ne divise pas le niveau, T p est 
un endomorphisme de l'espace. Lorsque p divise le niveau, on parle d'opérateur de 
Hecke spécial (et on écrit parfois U p , ce que l'on fera ici) ; dans ce cas, l'action sur les 
formes modulaires est un peu différente : elle est décrite par exemple par le lemme 1 
dans l'article de Li [5], et signifie intuitivement que l'opérateur fait "disparaître" une 
puissance de p du niveau, sans pouvoir faire disparaître complètement p. On verra plus 
loin une représentation graphique de la situation : on obtient une tour d'espaces de 
formes modulaires de dimensions finies, avec des inclusions croissantes avec le niveau, 
mais une dimension qui croît assez rapidement ; la projection est un moyen astucieux 
de contrer cette croissance. 

Expliquons rapidement la construction telle qu'elle est préséntée dans [7] : étant 
donnée une valeur propre a de U p agissant au rez-de-chaussée de la tour, on définit 
pour tout fi ^ le sous-espace caractéristique associé à (U p , a) : 

S fc Q W,X)= U Kcr(C/ p - ald)" 

permettant de définir une projection "K a ,^ via le choix d'un supplémentaire : 

n ImfC/p - ald)" 

la projection étant polynomiale en U p ; en faisant varier le niveau par des puissances 
de p, on obtient alors un diagramme avec projections -ïï a ._ '■ 



S k (Np^x)-^SZ(Np^ x ) 



S k (Np 2 lX )^^S%(Np 2 , x ) 



S k (N P ,x) 



S k (Np,x) 



dont la proposition 4.3 (de [7], toujours) affirme que les flèches verticales à droites sont 
des isomorphismes, d'où : n a .^ = (U, p M )~ 1 iT a ,iU! p . 



Ce dernier résultat est intéressant théoriquement, mais aussi pratiquement : on est 
ramené à appliquer successivement U p sur des espaces de plus en plus petits, puis pro- 
jeter au plus bas niveau, avant de finalement appliquer un inverse en petite dimension ; 
mieux, d'après l'article |2| de Coleman et Edixhoven, on sait pour k — 2 et condition- 
nellement à une conjecture de Tate pour k > 2 que l'opérateur U p est semi-simple au 
rez-de-chaussée, donc l'inverse en question est une simple multiplication par a -1 . On 
en voit un exemple concret en 14.21 

Remarquons d'ailleurs que cela signifie aussi que dans les étages supérieurs, tous les 
espaces caractéristiques associés à des valeurs propres non nulles de U p sont réduits aux 
espaces propres : la seule obstruction possible à la semi-simplicité de U p est la valeur 
propre zéro. Donnons un exemple concret : pour k = 2, N = 2,p = 5 et \ modulo 10 
vérifiant xCO = 1 ou x(7) = — 1, l'espace est réduit à zéro en niveau Np donc il n'y a 
pas de valeur propre non nulle, est un plan en niveau Np 2 sur lequel U p agit par zéro, 
et est de dimension 28 en niveau Np 3 ; U p est alors une matrice 28 x 28 de rang 2, dont 



ce qui montre que les conditions de semi-simplicité de [3] sont optimales. 
On prouve ici tout d'abord le résultat élémentaire suivant : 

Proposition 1 La projection caractéristique est un polynôme construit explicitement 
à partir du polynôme minimal de l'opérateur. 

(l'algorithme, utilisant l'algorithme d'Euclide étendu, est présenté de façon théorique 
en section |2~T| puis en détaillant plus les étapes algorithmiques en section l3~Tjl 
Les deux résultats d'indépendance suivants le rendent plus utile : 

Proposition 2 -Il existe toute une famille de polynômes convenables, qui donnent 
tous le même opérateur de projection (voir \2.ï\) . 
- L 'opérateur de projection ne dépend pas du choix du polynôme annulateur choisi 



Le premier point peut permettre de simplifier les calculs, au prix d'un polynôme 
de degré plus élevé (voir par exemple en section f3. 2 l une variante de l'algorithme avec 
moins de divisions) ; le second permet d'initier le calcul avec le polynôme minimal ou 
le polynôme caractéristique par exemple, ce qui donne un peu plus de liberté dans les 
calculs - mais ne les simplifie que de façon marginale ; le résultat le plus important 
pour les applications, qui permet des gains conséquents, est celui-ci : 

Théorème 1 La contruction de la projection canonique est équivariante ; plus précisément, 
étant donnés deux espaces vectoriels de dimension finie E et F, u G C{E) et v G C(F), 
j G C(E, F) et a £ C, tels que j o u — v o j , alors j o TiE,a{u) — ^F.a{v) ° j ■ 

(il est prouvé en 12 .4[) 

En effet, ce résultat a comme corollaire immédiat : 




(voir Wlfy . 



Corollaire 1 Dans l'égalité suivante, où p divise N, étant donné a € <D : 

S k (T 1 (N)) = ®S k (N,x) 
x 

on peut calculer les opérateurs de projection associés à (U p ,a) sur les sous-espaces de 
droite (plus petits), et obtenir par recollement l'opérateur de projection sur l'espace de 
gauche. 

Pour énoncer le second corollaire ci-dessous, il faut en dire plus sur les espaces 
de formes modulaires et expliquer un peu la théorie de Atkin-Lehner-Li : en plus des 
opérateurs de Hecke, on peut considérer les opérateurs Va, dont la définition est assez 
simple : pour d ^ 1, on a (Vdf){z) — f(d z) ou (Vdf)(q) — f(q d ), ce qui revient au 
même. En particulier, V\ est l'inclusion naturelle mentionnée précédemment. 

Si on fixe un niveau N ^ 1, on appelle forme ancienne de Sk(T\(N)) toute com- 
binaison linéaire de formes qui proviennent d'espaces de niveaux strictement inférieur 
via les opérateurs V - cela définit le sous espace Sk(Ti(N)) a . Le supplémentaire ortho- 
gonal de ce sous-espace (vis-à-vis du produit scalaire de Petersson) est appelé espace 
des nouvelles formes, noté Sk(Ti(N)) n . La théorie d'Atkin-Lehner-Li (développée par 
Atkin et Lehner dans pQ, puis par Li dans [5]) a pour conséquence que : 

S k (Ti(N) = e dM |JvK i 5 fe (r 1 (M)) n 

et cette décomposition est par exemple utilisée pour les calculs d'espaces modulaires 
par informatique dans SAGE, comme expliqué dans le livre de Stein [9]. 

On peut maintenant énoncer le second corollaire du théorème d'équivariance : 

Corollaire 2 Dans la décomposition précédente, si les seuls d qui interviennent réellement 
(pour lesquels les sous- espaces ne sont pas nuls) sont premiers avec p, alors on peut 
calculer les opérateurs de projections sur les sous-espaces de droite (de dimensions 
inférieures) et obtenir l'opérateur de projection sur l'espace de gauche par recollement. 

(l'exemple 14. Il illustre ce résultat) 

Je souhaite remercier Alexei Pantchichkine pour son soutien indéfectible, et William 
Stein pour m'avoir donné accès aux machines du réseau math. Washington. edu0. 

2 Généralités d'algèbre linéaire 

2.1 Définition élémentaire 

On se donne E un K-espace vectoriel de dimension finie et u un endomorphisme 
de E, admettant une valeur propre a. 

Le polynôme minimal de u sur E admet une factorisation de la forme (X — a) v Q 
où v ^ 1 et Q(a) / 0; d'après le théorème de Bézout, il existe alors un couple de 

2. Les calculs présentés dans cet article sont tous aisés et rapides sur une machine simple, mais les 
tâtonnements et expériences nécessaires à leur recherche ont parfois nécessité l'accès à une puissance de 
calcul plus conséquente ; les machines en question ont été financées par " National Science Foundation 
Grant No. DMS-0821725" . 



polynômes (A, B) tel que (X — a) p 'A + QB — 1, dont on déduit une décomposition de 
E en supplémentaires : E = Ker(u — a) v © Ker Q(u). 

L'espace Ker(it — a) v est appelé espace caractéristique de u associé à la valeur 
propre a ; c'est le plus grand espace sur lequel u — a est nilpotent. La décomposition 
en supplémentaires stables par u précédente permet alors de définir une projection sur 
cet espace caractéristique; on vérifie aisément que Q{u)B{u) donne une expression de 
l'opérateur associé. Cela établit déjà la proposition Q] 

Une remarque simple mais fondamentale est que le couple de polynômes (A, B) 
considéré n'est pas unique, mais que deux couples admissibles distincts diffèrent d'un 
multiple de (X — a) u Q - le polynôme annulateur de u. Donc l'opérateur Q{u)B(u) ne 
dépend pas du couple choisi ; ce qui prouve la première partie de la proposition p2J 

On notera ■nE,a{u) l'opérateur de projection caractéristique de u associé à a sur 
l'espace E. 

2.2 Stabilité par extension 

On se donne maintenant un K-espace vectoriel E de dimension finie, u un endo- 
morphisme de E ; et on suppose que F est un sous-espace vectoriel de E stable par u, 
et que a est une valeur propre de u sur F. 

On dispose donc de deux opérateurs de projection caractéristique, 7Te iQ (u) et TTF,a{u), 
que l'on souhaite comparer sur F. 

Les polynômes minimaux de u sur F et E sont de la forme respectivement (X— a) v Q 
et (X — a)^R, avec I ^ fi ^ v, Q diviseur de R et R(a) ^ 0. Écrivons R = QS pour 
fixer les notations. 

L'application du théorème de Bézout pour le polynôme minimal de u sur E fournit 
un couple (C, D) de polynômes tel que (X — a) v C + RD — 1 ; on a donc 7T£,a(w) = 
R(u)D(u). 

Mais cette même expression peut s'écrire (X — a)^(X — a) v ~^C + QSD = 1 ; le 
couple de polynômes ((X — a) v ~^C,SD) est donc admissible pour l'application du 
théorème de Bézout pour le polynôme minimal de u sur F ; on a donc 7Tf. q (u) = 
Q{u){SD)(u). 

Il reste à écrire Q{u){SD)(u) = (QS)(u)D(u) = R{u)D(u) pour conclure que l'on 
a bien 7i> jQ; (u) = 7TE, a (u)\F- 

Ce résultat est un lemme de base pour la preuve de la proposition[5]et du théorème[T] 

2.3 Changement de polynôme annulateur 

Une première variante importante de ce raisonnement est le cas où l'on ne considère 
u que sur un espace, mais avec deux polynômes annulateurs. Comme le polynôme mi- 
nimal est un diviseur commun à deux tels polynômes, il suffit de considérer la situation 
d'un polynôme annulateur quelconque par rapport au polynôme minimal : une com- 
patibilité dans ce cas montre alors que la compatibilité est générale. 

Maintenant, si compare avec la preuve précédente, on voit que la seule chose que 
l'on ait utilisée sur le polynôme minimal sur le sur-espace est qu'il était un multiple du 
polynôme minimal sur l'espace de base ; la même preuve s'applique donc : le calcul de 



la projection caractéristique ne dépend donc pas du polynôme annulateur considéré; 
cela établit la seconde partie de la proposition [2j 



2.4 Stabilité par morphisme 

Une seconde variante importante de ce raisonnement est celui où l'on considère 
deux espaces vectoriels de dimensions finies E et F, u un endomorphisme de E, v 
un endomorphisme de F et j une application linéaire injective de E dans F telle que 
j o u — v o j. Si a est une valeur propre de u sur E, c'est aussi une valeur propre de v 
sur F et on souhaite comparer 7r_E,a(w<) et 7Tp a (v) ; la situation peut se visualiser ainsi : 



Si P est un polynôme annulateur de u sur E, alors la relation jP(u) = P(v)j 
montre que P est un polynôme annulateur de v sur j(E) ; les mêmes raisonnements 
que précédemment montrent que le calcul de la projection caractéristique peut se faire 
au départ comme à l'arrivée : 



C'est le théorème d'équivariance[T] 

3 Pratique 

3.1 Calcul direct 

Il est assez facile de déterminer le polynôme minimal d'un endomorphisme donné 
sous forme matricielle via l'algorithme de Wiedemann, que l'on trouve détaillé en sec- 
tion 7.5.3 du livre de Stein [9]. 

Il est aisé de calculer un couple explicite dans le théorème de Bézout via l'algorithme 
d'Euclide étendu, tel qu'il est présenté par exemple en section 3.2 du livre de Cohen [2J. 

Le calcul d'un polynôme d'endomorphisme à partir d'un polynôme explicite et 
d'un endomorphisme donné par exemple sous forme d'une matrice, est aussi facile, via 
l'algorithme de Horner par exemple. 

Finalement, en combinant ces étapes, on obtient facilement une expression matri- 
cielle de la projection canonique ; la principale gêne dans les calculs est le passage au 
corps de nombres Q(a), dans lequel on réalise le calcul de divisions euclidiennes pour 
l'algorithme d'Euclide. 

Un cas particulier important est celui où l'on sait que a est une racine simple du 
polynôme annulateur considéré ; en effet, l'algorithme d'Euclide s'applique en une seule 
division euclidienne, que l'on peut calculer par la variante de l'algorithme de Horner 
qui traite la division par X — a ; elle ne fait donc intervenir que peu de produits et de 
sommes et une seule division. 

Un sous-cas intéressant du point précédent est celui de la tour modulaire au rez- 
de-chaussée de laquelle on sait que l'opérateur U p est semi-simple (sous couvert d'une 
conjecture de Tate pour k > 2), car alors on peut simplifier le polynôme annulateur 
choisi pour se ramener à une racine simple. 




j ° n E , a (u) = tt FiQ (w) oj 



3.2 Calcul sans division intermédiaire 



On sait que les divisions du calcul précédent sont incontournables en général. 

Cependant, les polynômes avec lesquels on travaille ont une forme assez particulière ; 
en utilisant cette spécificité, on peut contourner en partie cette difficulté de la façon 
suivante : étant donnée une factorisation d'un polynôme annulateur sous la forme 
(X - a) u P, avec P(a) = a ^ 0, on écrit P = a + (X - a)^Q avec ^ > 1. 

Pour e ^ 1, on a alors : 

Pf[(a 2 ^ + (-1Y(X - af-^Q*- 1 ) = a 2 ' - (X - af^Q 2 " 

de cette façon, si 2 e /z ^ v, ce qui se produit pour e assez petit, on a : 

pf[ (a 23 " + (-iy(X - af'^Q 23 ' 1 ) +(X- a )"(X - af^Q r = a 2 ' 
i=i 

expression qu'il suffit de diviser une fois par le scalaire a 2 pour écrire un couple ad- 
missible. On remplace le coût des divisions euclidiennes par un calcul avec un polynôme 
de degré plus élevé ; c'est un compromis. 

Remarquons aussi qu'une bibliothèque de calcul en théorie des nombres comme 
flintU dispose de routines optimisées pour les opérations de décalage (« Taylor 
shift >) et de mise au carré, qui peuvent être mises à profit pour les calculs précédents. 

3.3 Détermination des dimensions 

On a vu précédemment comment calculer très explicitement et complètement les 
opérateurs de projections caractéristiques, et on sait qu'alors la dimension de l'espace 
sur lequel on projette s'obtient par un simple calcul de trace. 

On peut cependant souhaiter connaître la dimension de l'espace sans avoir calculé 
la projection, par exemple parce que l'on n'a pas encore choisi quelle valeur propre 
utiliser. Dans ce cas, il suffit de calculer le polynôme caractéristique de l'opérateur : en 
effet, il suffira alors de le factoriser pour obtenir les ordres de multiplicité algébrique, 
donc les dimensions des espaces caractéristiques associés. 

4 Exemples concrets 

4.1 Calcul via la théorie d'Atkin-Lehner-Li 

On vérifie aisément (avec SAGE, qui permet de déterminer les espaces réduits à 
{0}) que l'on a la décomposition suivante : 

5 2 (ri(30)) = Fi5 2 (r 1 (15)) new © V 2 5 2 (r 1 (f5)) ne w © ^iS 2 (ri(30)) ncw 
(décomposition dans laquelle les espaces sont de dimensions respectives 9, 1, 1 et 7). 

3. http://www.flintlib.org/ 



Une base adaptée à cette décomposition est 



q-q 2 - ç 3 -q 4 + q 5 + q 6 + 3q 8 + q 9 + 0(q 10 ) 
q 2 - q 4 - q & - q 8 + 0(ç 10 ), 
q + 3q 8 -Uq 9 + 0(q 10 ), 
q 2 + 2q 8 -10g 9 + 0(q 10 ), 
q 3 + 2q 8 - 9q 9 + 0(q W ), 
q 4 + q 8 -5q 9 + 0(q 10 ), 
q 5 + q 8 -5q 9 + 0(q 10 ), 

q 6 -2q 9 + 0(q 10 ), 
q 7 + q 8 - 3q 9 + 0(q 10 ) 

La base calculée par SAGE directement pour le même espace est une base de Miller ; 
donc très simple : 

q + 0(q 10 ), q 2 + 0(q 10 ), q 3 + O(q 10 ),q 4 + O(q 10 ),q 5 + 0(q W ), 
q 6 + 0(q W ),q 7 + 0(q 10 ), q 8 + 0(q 10 ), q 9 + 0(q W ) 

La matrice de passage est donc : 



P = 



( 1 

1 






V o 



-1 1 

-1 







1 
1 





1 

-1 











1 
1 



1 \ 



-14 
-10 

-9 
-5 
-5 
-2 

-3 y 



L'endomorphisme T 3 agit par —1 sur les deux premiers sous-espaces, et par : 
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V o 





-3 
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o / 



sur le dernier ; son polynôme caractéristique est alors : X 7 + 3X 6 + 5X 5 + 7X 4 + X 3 — 
5X 2 — 3X — 9. Il a une racine simple réelle, 1, deux racines simples conjuguées dans 
(i et —i) et quatre racines simples conjuguées dans un corps de nombres de degré 
4 ; il s'agit donc aussi du polynôme minimal de l'endomorphisme. 

Choisissons de calculer la projection caractéristique associée à la valeur propre 1. 
Sur les deux droites, elle est donnée par zéro ; et sur le troisième espace, par le polynôme 



-L(X 6 +4X 5 + 9X 4 + 16X 3 + 17X 2 + 12X + 9) appliqué à la matrice précédente, soit : 



/ 



V 



1 1 

1 I 

1 I 

4 4 



1 1 

'4 4 



1 1 

"4 4 



■i i i 2 \ 

? ? ? ^ 

"i i i : 

4 4 4 1 







1 





.1 I I 1 / 

4 4 4 1 / 



Les résultats du paragraphe précédent ont été obtenus de la façon suivante dans 
SAGE : 



# Détermination de l'espace et de l'opérateur 
S30=ModularForms (Gammal (30) , 2) . cuspidal_subspace () 
S30n=S30 . new_subspace ( ) 

T=S30n.T(3) 

# Calcul du polynôme caractéristique 
R.<t>=PolynomialRing(QQ) 

P=T. charpoly(t) 

# Calcul du polynôme à utiliser pour le calcul de la projection 
Q=P.quo_rem(t-l) [0] 

Pproj=xgcd(Q,t-l) [1]*Q 

# Obtention de la matrice de la projection 
Pproj (T.matrixO) 

La matrice de la projection canonique sur l'espace complet est alors, dans la base 
adaptée à la décomposition donnée ci-dessus : 
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'4 4 



"1 î î 9 

"! ! ! i 

"î î î ; 

4 4 4 L 



.1 I I 1 

4 4 4 ^ 


111 

'4 4 4 



1 / 



En calculant directement avec l'opérateur sur l'espace total, il faut appliquer sa 
matrice au polynôme ^ X s + 6X 7 + 18X 6 + 38X 5 + 58X 4 + 62X 3 + 50X 2 + 'S0X + 9, 



d'où : 



3 «v "V 3 3 _ 2 7 3 «V 

I _I I I _I _I _9 _I I 



| _| | | _| _| _2 -| | 




matrice qui est équivalente à celle trouvée précédemment, à P près. 



4.2 Calcul avec une tour modulaire 

Cet exemple est un cas où les calculs sont encore possibles à tous les niveaux, mais 
illustre le gain du passage à une dimension inférieure. 

On choisit un niveau de base 30 et le nombre premier 3, et un caractère de Dirichlct 
trivial, de sorte que la situation est : 



(dimension 11) <S 2 (90,1) ^S^ 1 (9Q,1) 



(dimension 3) 



u 3 



u 3 



5 2 (30,l)^ r7t 5 2 - 1 (30,l) 

Us Us 



On constate simplement que U3 agit par simple multiplication par — 1 sur les espaces 
caractéristiques, qui sont des plans ; la tour permet donc de calculer la projection en 
haut via un calcul en bas, en dimension inférieure. 

Ici, la projection au rez-de-chaussée se calcule avec le polynôme — jX 2 — + | 
et est donnée par la matrice : 

1 1 

2 2 

1 

1 1 

2 2 

A l'étage au dessus, elle est donnée par le polynôme — ^-X w — ^X 9 + ^-X 8 , et la 




matrice 
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V 

En choisissant une forme modulaire / dans l'espace ^2(90, 1), on vérifie aisément 
l'égalité 7T-i,x(/) = — 7T-i.oC^3(/) ! l a session ci-après utilise les deux fonctions acces- 
soires données en section [5] : 

sage: # Chargement des fonctions d'aide 

sage : load( ' /home/ jpuydt /Recherche/abaisse . sage ' ) 

sage: load( ' /home/jpuydt/Recherche/polyproj .sage') 

sage: # On travaille avec des polynômes, pas des expressions 

sage: R. <X>=PolynomialRing(QQ) 

sage : # Données du calcul 

sage: N=10 

sage : p=3 

sage : N0=N*p 

sage : # Choix du caractère 

sage: chi=DirichletGroup(NO) [0] 

sage: # Calcul de la projection en bas 

sage: Eb=ModularForms (chi) . cuspidal_subspace () 

sage: Ub=Eb.T(p) 

sage : polb=polyproj (Ub . charpoly (X) ,-1,2) 

sage: # Calcul de la projection en haut 

sage: Eh=ModularForms (chi . extend(N0*p) ) . cuspidal_subspace () 

sage: Uh=Eh.T(p) 

sage : polh=polyproj (Uh. charpoly (X) ,-1,2) 

sage : # Comparaison des deux 

sage: f=Eh.gens() [2] 

sage: g=abaisse(Eb, Uh(f)) 

sage: ordre=Eh . sturm_bound () 

sage: polb(Ub) (g) .qexp(ordre) + polh(Uh) (f ) .qexp(ordre) 
0(q~37) 

sage: # On obtient bien zéro! 



5 Code SAGE 



5.1 Fonction polyproj 

def polyproj (poly, racine, ordre): 

"""Reçoit en argument un polynôme annulateur d'un endomorphisme , 
une de ses racines et son ordre et retourne un polynôme permettant 
de calculer la projection caractéristique sur la racine 
correspondante " " " 
var=poly. parent () .gen() 
pow= (var-racine) **ordre 
Q=poly . quo_rem(pow) [0] 
return Q*xgcd(Q ,pow) [1] 

5.2 Fonction abaisse 

def abaisse (espace, forme): 

"""Reçoit en argument un espace de formes modulaires et une forme 
modulaire a priori dans un espace de niveau plus élevé, mais dont 
on sait qu'elle devrait être plus bas, et retourne la 'même' forme 
en bas. 

coef f s=espace . f ind_in_space (forme) 

return espace . linear_combination_of _basis (coef f s) 
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